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This  report  relates  the  characteristics  of  the  Holloman  High  Speed  Test 
Track's  simulated  rainfleld  to  those  of  a  natural  rainfleld  as  defined  by 
mathematical  models  propogated  by  other  experimenters.  Data  collected  for 
heavy  rain  rates  (25-30  In/hour)  generated  for  Project  Deluge  In  a  controlled 
environment  are  analyzed  and  a  tentative  mathematical  model  relating  the  key 
parameters  Is  provided.  This  report  specifies  the  additional  efforts  required 
to  finalize  the  model. 
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FOREWORD 


The  purpose  of  rain  erosion  testing  on  the  High  Speed  Test  Track  (HSTT)  at 
Holloman  Air  Force  Base  is  to  study  the  erosive  effects  of  extended  supersonic 
or  hypersonic  flight  through  rain  clouds  on  material  samples  and  components 
of  aerospace  systems. 

The  mechanism  used  to  generate  the  simulated  rain  environment  is  a  series 
of  adjustable  spray  nozzles  spaced  along  the  side  of  the  track  adjusted  to 
discharge  over  the  flight  path.  The  rain  rate  is  controlled  by  the  water 
pressure  and  drop  size  distribution  is  controlled  by  a  combination  of  nozzle 
adjustments  and  pressure. 

For  the  selection  of  materials  to  be  effective,  it  is  essential  that  the 
component  designer  understand  the  relationship  of  the  simulated  rainfield 
characteristics  and  those  of  natural  rainfields.  This  report  is  the  first 
attempt  to  quantify  these  relationships,  and  provide  the  designer  with  a 
tool  to  assess  the  differences. 
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PROJECT  DELUGE:  A  PROPOSED  MATHEMATICAL  MODEL  FOR  THE  HOLLOMAN  AFB  TEST  TRACK 
SIMULATED  RAINFIELD,  VOLUME  II 

1.0  BACKGROUND 

1.1  During  the  performance  of  Project  Deluge  it  became  apparent  that  a 
requirement  existed  for  a  means  to  define  the  Test  Track's  simulated  rainfield 
capabilities  to  achieve  the  parameters  most  often  used  by  prospective  customers 
to  identify  their  requirements,  i.e.  accumulated  rain  rate  (ARR)  and  mass  median 
diameter  (MMD). 

1.2  Prior  research  of  natural  rain  has  demonstrated  that  the  drop  size 
distribution  of  raindrops  in  a  rainfield  can  be  defined  by  an  exponential 
expression  of  the  form  (Ref  1) 

Nj  =  No  EXP  {-  J{,0) 

where: 

Nj  =  number  of  drops  of  unit  size  range  per  unit  volume  (N/m  ) 

N  =  Value  of  N.  when  D=  0 
0  d 

D  =  Drop  diamter  (mm) 

_A,  =  Variable  which  is  dependent  upon  rain  rate  or  intensity  of  precipitation 
The  relationship  of  A  with  rain  rate  is  defined  by  the  expression: 
aR^  (mm'^ ) 

where:  a  and  b  are  constants 

R  =  rain  rate  (mm/hr) 

1.3  The  following  methods  were  used  by  past  experimenters  to  derive  the 
relationships  noted  above: 

a.  Rain  rate  was  determined  by  rain  gauge  measurements 

b.  Drop  size  distributions  were  determined  by  using  flour  beds  to 
catch  rain  and  calculating  the  size  from  sieves  or, 

c.  Measuring  the  impressions  of  drops  on  blotting  paper 

1.4  A  parallel  methodology  has  been  employed  to  formulate  the  equations 
for  the  Test  Track  simulated  rainfield.  The  major  difference  is  the  use  of  a 
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mass  spectrometer  to  count  and  catalog  the  raindrops  by  unit  size.  The 
substitution  of  the  spectrometer  for  flour/blotting  paper  win 
enhance  the  accuracy  In  determining  the  statistics  for  the  drop  size 
distribution  and  the  resulting  equations. 

2.0  DATA  BASE  The  data  used  In  this  report  was  accumulated  as  part  of 
Project  Deluge  during  November  1981.  The  mass  spectrometer  used  to  measure 
and  count  the  drop  sizes  was  the  Illinois  Institute  of  Technology  Research 
Institute  (IITRI)  Counter  specifically  designed  and  constructed  for  use  In 
the  Test  Track's  simulated  rainfleld.  The  use  of  this  counter  is  discussed 
In  Volume  1,  along  with  the  test  conditions,  raw  data  formating  and  data 
output  from  the  counter.  (Reference  AD-TR-83-42,  The  Development  of  a  Holloman 
High  Speed  Test  Track  Heavy  Rainfleld) 

2.1  The  basic  parameters  measured  during  the  experiment  were: 

a.  Rain  rate  -  using  standard  rain  accumulation  type  gauges  (ARR)(1n/hr) . 

b.  Drop  size  and  number  of  each  -  IITRI  Counter.  Drop  size  interval 
width  of  0.6mm,  centered  around  diameters  of  0.6,  1.2,  1.8,  2.4,  3.0,  3.6,  4.2, 
and  4.8mm. 

2.2  Data  parameters  calculated  and  provided  by  Project  Deluge  are: 

a.  Number  of  drops  per  unit  size  per  cubic  meter  (N/m  ) 

b.  Mass  Median  Diameter  (MMD  or  Do),  (mm). 

c.  Liquid  Water  Content  (LWC),  (grams/m  ). 

d.  Equivalent  Rain  Rate  (ERR),  (In/hr). 

3 . 0  DATA  ANALYSIS  AND  EVALUATION 

3.1  The  data  of  paragraph  2.0  above  was  subjected  to  regression  analysis  to 
determine  the  best  fit  relationships  between  the  different  variables.  The 
computer  program  selected  for  this  function  provides  the  capability  of  evaluating 
seven  regression  models  simultaneously  for  each  set  of  data.  The  resulting 
equations  providing  statistically  significant  relationships  are  provided  below 
with  the  result  of  previous  experimenters. 

3.2  The  matrix  format  provides  a  means  of  comparing  previous  experimenters 
results  with  natural  rainfall  and  the  results  obtained  with  Project  Deluge. 
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Type  of 

Source  Precipitation 


3.3  Other  Established  Relationship: 
Dyer  ( Ref  1 ) 

Project  Deluge 


Do  =  3.67  A 
Do  =  3.6779  A 


-1 .00 
-0.996 


Coefficient  of  Determination  =  0.875 


3.3.1  Assessment  of  Statistical  Parameters: 

Coefficient  of  Correlation  -  The  hypothesis 


P=  0 
P  >  0 

was  tested  using  the  student's  t  test,  where  the  test  value  is  defined  by: 


t'-iV 


N-2 


with  acceptance  criterion 


■*- 

t 


and  with  rejection  criterion 


★ 


t  ^ 


*N-2 


By  manipulating  the  equation  for  t  ,  it  can  be  shown  that  for  a  sample  size  of  81, 
the  following  minimum  values  of  r  is  required  for  rejection  of  H^:  at  the  levels 
of  significant  )  Indicated. 


or 

0.10 

0.05 

C.02 

0.01 

O.l'^PO  j 

0.21 53 

0.2532 

0.2783 

-1  -3 

(a)  Nq  (mm  m  )  -  Defined  as  the  number  of  drops  per  unit  size 
per  unit  volume  when  droplet  diameter  is  zero.  Mathematically, 

Limit  Un  = 

D  0 


For  the  Marshall-Palmer  model  this  is  defined  to  be  a  constant  value  of  8000. 
The  best  relationship  derived  fromthe  simulated  rainfield  is: 

=  R  0.117  +  2.288X10’^  R 

This  equation  implies  that  as  rain  rate  (R)  goes  from  zero  to  infinity, 
goes  from  zero  to  43706.3.  Mathematically  given  by: 

(1 )  Lim  =  Lim  R 

R--?0  R-0  Q  1^7  ^  2.288X  10'®  R 

0 

0.117  +  2.288X  10-5  (0) 


Lim  Nq  =  0 

R"© 


{i)  Lim  =  Lim 

R 

0.117  +  2.288X  10'^  R 

Lim 

1 

0.117  ^ 

R 

2.288X10"® 

1 

0  +  2.288X10"® 

L  im  N 

0 

43706.3 

R 

Conclusion:  There  is  significant  difference  in  both  the  value  and 
characteristics  of  between  the  two  models. 
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(b)  LAMDA,  (  A  »  A  variable  which  Is  dependent  tipoo  rain  rate 

Generally  stated  mathematically  as: 

A  =  a  R  ^ 

For  the  Marshall -Palmer  model  the  values  of  a  and  b  are  4.1  and  -0.21,  respectively. 
The  calculated  values  of  a  and  b  for  the  simulated  rainfield  are  3.346  and  -0.125 
with  standard  error  of  estimates  of  1.146  and  0.021,  respectively.  Using  these 
values  the  Indicated  confidence  limits  were  calculated  to  be: 


(!-.<  )  % 

a  =  3.346 

b  =  -0.125 

99.00 

0.389  Sa  <6.302 

-0.1794  Sb  S  -0.0711 

95.00 

1.100  Sa  S  5.592 

-0.1664S  b  S  -0.0841 

90.00 

1 .461  Sa  <  5.231 

-0.1598  bS  -0.0907 

z  c 

0.658 

4.04 

Conclusions:  The  mathematical  relationship  betweenA  ^nd  R  is  confirmed  to  be 
a  power  function.  While  there  is  no  statistically  significant  difference  in  the 
scaling  parameter  (a)  values, there  isa  very  significant  statistical  difference 
in  the  shaping  parameter's (b)  values. 

(c)  Mass  Median  Diameter--Do(mm) ,  defined  as  the  droplet  -  diameter  which 
divides  the  Liquid  Water  Content  (LWC)  into  two  parts.  Generally  defined 
mathematically  by  the  equation: 

Dq  =  a, 

For  the  Marshall -Palmer  model,  the  values  of  a^  and  b  are  0.90  and  0.21, 
respectively.  The  calculated  values  of  a^  and  b  (simulated  rainfield)  are 
0,916  and  0.154  with  standard  error  of  estimates  of  1.143  and  0.0206,  respectively. 
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Again  the  indicated  confidence  limits  are: 


(l-« )% 

a^  =  0.916 

b  =  0.154 

99.0 

-2.033  S  a^  S  3.8659 

0.100  S  b  S  0207 

95.0 

-1.324  S  a^  ^  3.157 

0.1132  SbS  0.194 

90.0 

-0.9643  S  a^  S  2.797 

Z 

c 

0.0143 

2.734 

Conclusions:  Again,  the  basic  relationship  between  the  variables  is  confirmed. 
Also,  there  is  no  statistically  significant  difference  in  the  values  of  a-j ,  but 
there  is  a  statistically  significant  difference  in  the  b  values.  However,  note 
that  the  critical  Z  value  is  less  than  3.0,  i.e.  the  Marshall-Palmer's  specified 
b  value  is  less  than  3  (f  from  the  value  determined  for  the  simulated  rainfield. 

(d)  Overall  conclusions:  The  expected  relationships  for  both  JC  and 
as  a  function  of  rain  rate  were  verified.  The  significant  differences  in  the 
regression  coefficients  and  the  form  of  N^,  simply  implies  that  the  simulated 
rainfield  as  configured  does  not  duplicate  the  droplet  diameter  distribution  of 
an  Ottawa  summer  rain. 

4.?  Vs  Jones  (1956)  Model:  This  model  was  postulated  from  observations 
of  thunderstorms.  Definitions  of  the  parameters  are  the  same  as  previously 
provided.  Comparing  the  simulated  rainfields  to  Jones*  provides  the 
following  information: 

1  OP 

(a)  Nq  (mm  m  ).  Jones  values  is  =  3.89R  vs 
-  R  [  0.117  +  2.288X10'^lj 

For  both  situations  =  f{R),  with  only  the  form  of  the  relationship 
differing.  If  the  simulated  rainfield  data  is  fitted  to  a  power  function  to  place 


it  into  the  same  form  as  the  Jones'  observation,  the  fo  equdi  on  .  e:.. 

0.3959 

-  1925. 73R 


The  confidence  limits  for  the  regression  coefficients  for  the  the  power  function  are 


a^  =  1925.73 

b  =  0.3959 

99.0 

1816.35  S  a  S  2035.11 

95.0 

1842.64  S  a  S  2008.82 

1  .  - 

0.2489  SbS  0.5428  i 

90.0 

185.99  S  a  S  1  995.47 

0.2725  S  b£  0.5192  ' 

36.2 

. -  — ^  ' - - 1 

g.3  ( 

i 

1 _ : 

Conclusions:  Both  of  the  regression  coefficients  are  statistically  significantly 
different  from  those  identified  by  Jones.  Also,  although  both  situations 
identify  as  a  f(R),  the  forms  of  the  relations  are  different. 

(b)  A  (mm‘^)  .  The  Jones'  determined  relationship  is  A=  2.5R'^‘^^ 

vs  Track's  3.346R  jhe  table  of  confidence  limits  provided  in  4.1(b)  above 

is  applicable  with  the  same  conclusions.  The  critical  Z  scores  for  a  and  b, 
are  0.74  and  3.6,  respectively. 

(c)  Dp (mm).  Jones  1.48R°-°®  vs  Track  0.91 

The  Table  of  Confidence  Limits  provided  in  4.1(c)  above  is  applicable  and  the 
conclusions  are  the  same.  The  critical  Z  scores  are  0.493  and  b.05  for  a  and  b 
respectively. 


9 


» 


(d)  Liquid  Water  Content.  HwL  (q/m^) .  Defines  the  liquid  water  content 
of  the  rainfield  per  unit  volume.  Using  the  Jones'  equations  Oyer  defined  the 
following  relationship  for  LWC  as  a  function  of  rain  rate: 

(1)  J]lWC  =  0.052R°'^^ 

Independently,  Atlas, using  the  Marshall-Palmer  equations  derived  this  relationship: 

(2)  ^LWC  =  0.072R°'®® 

The  relationship  between  these  two  variables  observed  during  DELUGE  is: 

(3)  ^LWC  =  0.07889R°’®^®  (Regressing  R  on  ^LWC) 

"he  confidence  limits  for  the  derived  Deluge  values  are: 


(l-«  )% 


a  =  0.07889 


h  =  0.868 


99.0 


-4.458  Sa  ^  4.616 


0.644  <  b  SI. 0926 


95.0 

90.0 


-3.368  Sa  £  3.526 


-2.814  Sa  S  2.972 


0.698  S  bS  1 .0387 


0.725  S  bS  1 .0114 


2  Critical 


0.004 

0.015 


0.138 
1 .169 


(Atlas) 

(Dyer) 


Cone  I  I.-'"  )n:  There  are  no  statistically  significant  differences  between  the 
coeff'f  ienti  for  the  relationships  described  by  the  equation: 

^LWC  -  aR*^’ 

As  i.  by  '  cher  Atlas  or  Dyer  and  those  experimentally  determined  by  Track. 

.c  Overall  Conclusions:  The  expected  relationships  for  all  of  the 
i  cr  jf. ' '‘ir  .  were  verified.  The  differences  in  the  coefficients  are  most  probably 
attributable  to  differences  in  the  rain  rates  observed.  Deluge  rain  rates  of 
?9  inc^e'.  .'  e.  v/ould  be  observed  only  on  rare  occasion  in  nature. 
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4.3  Vs  Joss  et  a1  (1 970) .  This  model  is  Od:,e.  ti.jiu ur  ;. 
observations  and  is  very  similar  to  the  Marshal  1 a i.:ier  model.  The  differences 
between  the  two  models  (Joss  vs  Marshall)  are  the  values  cited  for 

(both  constants)  and  values  for  scaling  parameters  for  both -A.  and  as  f(R). 

The  comments  and  conclusions  drawn  for  the  comparison  of  Marshall-Palmer  vs 
Track  are  applicable. 

4.4  Vs  Dyer.  In  her  investigation  of  the  other  experimenters*  work 
Dyer  developed  the  following  relationship  between  ^  and  i.e.: 

A  =3,67 

or  Dg  =3.67A'^ 

This  relationship  appears  to  be  consistent  regardless  of  either  the 
characteristics  of  the  rainfall,  i .e., thunderstorm,  summer  shower,  etc.  or 
the  investigator.  This  relationship  also  holds  for  the  simulated  rainfield. 

The  values  derived  via  regressions  are: 

Ojj  =  3.6779  A 

It  can  be  readily  shown  that  no  statistical  difference  exists  between  the  equation 

4.5  Proposed  Model  for  Mini-Rainfield,  The  following  equations  are 
tentatively  suggested  for  use  in  describing  the  mini -rainfield  when  it  is 
configured  as  specified  in  the  Project  Deluge  report: 

General  Equation:  EXP  {  '  A  D) 

=  R  [0.117  +  2.288X10‘^r] 

D  =  0.91 
0 

Dq  =  3.6779A'' 

Rain  rate  (R)  is  the  measured  value  for  rain  rate  in  mm/hr. 

NOTE:  (1)  An  alternate  equation  for  is:  1 925.7R®’^^^^ 

(2)  J^LWC  may  be  estimated  using  this  equation: 

J^LWC  =  0.07889R°'®^  (g/m^) 


11 


RECOMMENDATIONS  FOR  ADDITIONAL  TESTS: 


In  some  instances  where  statistically  significant  differences  exist,  most  can 
De  attributed  to  the  fact  that  the  rain  rates  produced  for  Deluge,  i.e.  25-30 
inches  per  hour  are  almost  never  found  in  nature. 

While  it  has  been  shown  by  the  experimenters  referenced, that  the  parameters 
of  interest  are  dependent  upon  rain  rate  for  the  rain  rate  intensities  observed, 
the  sensitivity  of  the  indicated  parametric  relationships  to  changes  in  rain 
rate  has  not  been  identified. 

The  simulated  rainrield  at  Test  Track  has  the  capability  to  produce  a 
w''de  spectrum  of  rain  rate  intensities  to  allow  such  a  sensitivity  analysis. 

Tne  following  additional  tests  are  planned: 

{1}  To  extend  the  applicability  of  the  proposed  model  from  a  "specific 
condition"  of  the  mini-rainfield  to  a  general  statement  of  the  mini-rainfield 
capabilities,  additional  measurements  are  required  with  the  water  pressure  values 
set  at  vatues  incrementally  spaced  between  the  maximum  and  minimum  allowed. 

(2)  Because  the  Knoll enberg  Probe  has  the  inherent  capability  to  provide 
better  resolution  of  droplets  diameter,  that  it  be  used  to  acquire  future  data 
or  the  data  for  each  set-up  be  acquired  with  both  the  Probe  and  the  Counter. 

(3;  All  of  the  above  tests  will  be  repeated  on  the  Track  so  that  the 

effects  of  out-of-door  conditions  can  be  incorporated  Into  the  final  model  equations. 

The  selection  of  the  final  model  for  the  Test  Track  simulated  rainfield 
facility  will  be  decayed  until  after  a  comparison  is  made  between  the  two 
counters  and  the  abov^'  planned  tests  are  completed. 

T  ^>0rfJM^NTS- 

1.  Extremes  of  Hydrometers  at  Altitude  for  MIL-STD-210B,  Supplement  Drop  Size 
Distrloutiuns  ,  ATCrjl  -TR-71-0008,  9  Jan  73, 

2.  Smrtcr  tontr i hution  (Optical  Extinction  by  Rainfall),  Atlas,  David,  Journal 
of  Metecroicqy,  Vol  10,  Occ  1953. 
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